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ABSTRACT 


An. initial study has been completed to begin development of a flight 
experiment to measure spacecraft control/flexible structure interaction. 

The work reported consists of two phases: identification of appropriate 

structural parameters which can be associated with flexibility phenomena, 
and suggestions for the development of an experiment for a satellite con- 
figuration typical of near future vehicles which are sensitive to such 
effects. Recommendations are made with respect to the type of data to be 
collected and instrumentation associated with these data. The approach 
consists of developing the equations of motion for a vehicle possessing a 
flexible solar array, then linearizing about some nominal motion of the 
craft. A set of solutions are assumed for array deflection using a con- 
tinuous normal mode method and Important parameters are exposed. Inflight 
and ground based measurements are distinguished. Interrelationships between 
these parameters , measurement techniques , and input requirements are dis- 
cussed which assure minimization of special vehicle maneuvers and opti- 
mization of data to be obtained during the normal flight sequence. Limited 
consideration is given to flight data retrieval and processing techniques 
as correlated with the requirements imposed by the measurement system. 
Results Indicate that inf, light measurement of the bending and torsional 
mode shapes and respective frequencies, and damping ratios is necessary. 
Other parameters may be measured from design data. Areas of proposed 
investigation include examination of specific array configurations, com- 
puter simulation of array dynamics, interaction properties of the attitude 
control system, and design of the. measurement system. 
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mode dependent Integral expressions defined in 
Appendix B 
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frame 

k 

disturbance values for angular velocity of vehicle 

position vector (1 = 1,.. .6) 

time 

transformation matrix between {b} and (a) 
transformation matrix between {b} and {i} 

transformation matrix components 

force distribution on array differential element 

force distribution on array 

force magnitude on appendage element 

component or F s in a 3 direction 

f f F 8 <J, i (r 3 )dr3dr 4 
r 3 4 


R 0 7 4-121 


x 

p 

M 

N 

R 

R' 


dT 


a 





T 


1 


(Qi 0 a» 


CA(X))j 

U 


2 s 

6 1. 

X 

+i 

*i 


vili 


angular momentum vector 
moment of inertia of array element 
vehicle moment of inertia (p » 1,2,3) 
vehicle total mass 

degree of highest harmonic mode retained in analysis 

/ dr 4 
r 4 4 

/ r.dr, 
r 4 4 
4 

torque acting on array differential element 

s 

component in direction of T 

S I f ,♦ i (r 3 )dr 4 dr 3 
r 3 4 

angular velocity amplitude, frequency dependent 


fch 

bending amplitude of j mode 

th 

torsional amplitude of J mode 
mass ratio 

array torsional dependent variable 

array element rotation vector 

Kronecker delta 

frequency parameter 

bending mode shape of i^ 1 mode 

til 

torsional mode shape of i mode 


/ 


1 


R 0 74-121 



X. INTRODUCTION 

1. PURPOSE OF THE INVESTIGATION 

Adverse effects of flexibility on satellite attitude motion were \ 

i 

manifested early in the U.S. space program and have remained Important - 

factors in spacecraft designs. Mathematical analysis supported by ground ' 

experimentation have been relatively reliable for the present generation 
of satellites. However, future and presently developing spacecraft are 
progressing toward larger and highly flexible antennas and solar panels. 

In such cases ground experiments are impractical because of size and 

environmental restrictions, } 

The overwhelming Importance of the flexibility problem coupled with 

I 

possible inadequacies of present analytical approximations and ground test- 
ing procedures indicate the introduction of a program and apparatus 
capable of yielding inflight data on appendage dynamics and their effects 
on craft control. Only by implementation of a monitoring system on these 
new satellite configurations can continued progress be made on the control 
of these effects. In addition, further refinements on design procedures 
can be derived from this data. 

It i9 the purpose of this investigation to identify the structural 
and control system parameters which can be associated with the flexibility- 
control interaction phenomenon, and develop an experiment for a satellite 
configuration typical of near future vehicles . 

2. CURRENT AND RECENT FLEXIBILITY PROGRAMS 

Many prc /bed vehicles have appendages which are very large and | 

highly flexible. Examples are the radiotelescope satellite (Ref. 1,2), 
composed of flexible aluminum ribbons having a diameter of fifteen hundred 
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meters. Other equally Illustrative examples are abundant, i.e., the pro- 
posed direct T. V. broadcast satellite (Ref. 1) which possesses solar 
arrays with an area of one thousand square meters. The size of these 
vehicles indicates the importance of accurately anticipating flexibility 
interactions. 

Programs designated for flexibility experiments have generally been 
among the first to suffer frov\ cost cuts, weight decreases, and other 
considerations leading to system abolishment. Minor flexibility control 
interaction tests were performed in conjunction with the Gemini, Apollo 
programs (Ref. 3,4), and a major attempt at flexibility measurement was 
planned in conjunction with the Canadian Communications Technology satel- 
lite (Ref. 5,6). However, this effort has been stripped to an optical 
system which will measure only static deformations. The Flexible Rolled 
Up Solar Array (Ref. 7) has been flown successfully and is equipped with 
six accelerometers and three strain gages and temperature sensors for the 
detection of array dynamic response. Interaction effects were not observ- 
able and peak accelerations have been only of the order of a few milli-g's. 

Several significant analytical studies have been undertaken to deter- 
mine appendage dynamics and control interaction phenomenon. Studies have 
been completed which Include the contribution of large flexible appendages 
to control difficulties arising from continuous and flexibly connected 
bodies. One of the more promising analytical methods uses hybrid coordi- 
nates (Ref. 8,9). This method formulates the flexible space vehicle con- 
trol as a combination of discrete and modal coordinates. The analysis 
involves three distinct steps. First, preliminary design is based on root 
locus plots for single axis response of linearized systems. Second, modi- 
fications are made as requii'ed by eigenvalue analysis of the coupled 


.linear systems. Third, nonlinear differential equation simulation is 
accomplished using numerical integration. The third step is used for 
design confirmation. 

In particular cases, several other methods of analysis are adequate 
depending on vehicle geometry and appendage characteristics. These are 
the energy sink method, usually applicable to the compact, near rigid 
body type of vehicle; the discrete parameter method, which is best suited 
to a satellite with compact bodies and flexibly attached appendages; and 
the modal coordinate method, which is adaptable to highly flexible craft 
with larger overall dimensions, 

X! . FLEXIBLE STRUCTURE 

1. DEFINITION 

An adequate definition for flexibility may be evaluated only with 
regard to the context in which it is to be used. For applications to 
satellite flexibility control interaction, it is not necessary to define 
flexibility in the strictest physical sense but rather to define those 
aspects of flexibility necessary in predicting important coordinate values 
and parameters of the interaction. It is also necessary that the chosen 
mathematical model generate equations of motion which are computationally 
treatable. 

Two realms of flexibility-control interaction immediately distinguish 
themselves as important. First, attention should be directed toward deter 
mining attitude control sensitivity to the elastic mode frequencies of 
appendages and flexibly connected subsystems of rigid bodies. Secondly, 
sensitivity of attitude control sensors to changes in vehicle Inertial 
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properties arising from flexibilities must be investigated. Restricting 
the definition of "critical flexibility" to those considerations critical 
to the interaction problem leads to the following definition! 

Critical flexibility with respect to the attitude control- 
structural dynamic Interaction problem are those modes of 
flexibility which create elastic deformations of continuous 
and discrete masses with modal frequencies within the fre- 
quency sensitivity range of the control system, and also 
those flexibility states which cause redistribution of 
vehicle Inertia to an extent that control system activa- 
tion is generated. 

HI. DEVELOPMENT OF EQUATIONS OF MOTION 

1. ASSUMPTIONS AND MODEL CONFIGURATION 

The spacecraft shall be modelled using a relatively simple configura- 
tion. It is found that this greatly reduces mathematical complexity of 
the model without significantly compromising the utility of resulting 
equations. The configuration chosen consists of a. rigid central body with 
one continuous flexible solar array. The second array may be considered 
stiff. The appendage axis in its undeformed position Is assumed to act 

through the vehicle center of mass. It is assumed that no external forces 

or torques are applied to the array and that the control system will react 

with a pure torque only. The basic configuration of the craft, earth, sun 

system is shown in Fig. 1. 

Three coordinate systems are used in the formulation of the problem. 

An inertial reference {i} is fixed at the earth's center. A second system 
{&} is fixed to the vehicle center of mass and is allowed to rotate with 
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the rigid central, body. The third coordinate system {a} is fixed in the 
solar array and allowed to rotate with the array as it maintains proper 
solar orientation. The origin of the {a} system will coincide to the 
origin of the {£} system when the array is undeflected, as will the 
and bj axes. For further simplification, the array is assumed flat and 
of homogeneous construction with no mass discontinuities. This assumption 
greatly reduces the analytical complexity while retaining the presence 
of important parameters. Array position vectors are shown in Fig. 2. 


2. ARRAY EQUATIONS OF MOTION 

If an element of the array is considered as shown In Fig. 2, the 
position vector is the sum of those vectors which may be referred to the 
pruviously defined reference frames. The equation of translational motion 
for the element may be derived from Newton's second law. 



•+ 

» t r .dm 
i-i 1 


a) 


The frames of reference may be related through the use of appropriate 
transformation matrices. 

{a} « C 1 ' {b} 

and 

fb} = C 2 {i} 

The motion of the {a} frame with respect to the {£} frame will con- 
sist of a rotational motion known in time and a translational motion which 
Is a function of the array deformation. For small time intervals, the 
motion of {t>} with respect to {£} may be expressed in terms of the rigid 
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body angular velocity variables q^. Transformation components of 0 and 

2 12 
C shall be designated and . 

Disturbance quantities shall be expressed as follows. Each component 

of a vector 1, shall consist of a nominal steady state value (e.g. A^) 

plus a disturbance quantity (e.g. A^q) . In reference frame {a}, for 

example, vector will appear as, 

K ° Ca l a 2 a 3 ] 

where the superscript specifies the vector component. Position vectors 
are specified in Appendix A. 

Using Eq. (1), a matrix relation extracted from work by Likin's 
Hef. 8), gives for the array element; 

F S " m 8 {(^[C 2 ^ + r 2 + 2q x r 2 + t^q^ + q^] 

^ /-w, 6 

- ( S r 1 )‘"(C 1 q 1 ) + [ (C^) (G 1 q 1 ) + 2(C i q 1 )q 2 -I- q 2 q 2 I< * r ±) 
i=3 

+ 2[qf 2 + (C 1 q 1 )]<r 5 + r 6 ) s + (r 5 + r fi ) s } (2) 

where the coordinate frame is ia). 

The vector r 2 exists only due to the array deformation. Therefore, 

-h . 4- 

it may be expressed in terms of the deflection variables r^ and r^. 
Defining the terms. 


(Ai + A 1Q ) 

(A 2 * A 2 ,) 

(Ai + A 10 ) 


K o t ~ & r 
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m 

M 


R 


dr. 


(3a) 

(3b) 


-I 


R - ' r 4 dr 4 

r 4 


(3c) 


allows expression of as > 


r 9 » -]h 


C Z [r 5 R + JR' ]di ; 3 


(4a) 



-P 


C 2 [r 5 R + oR' ]dr 3 
3 


(4b) 



| C 2 [r 5 R + ciR f ]dr 3 
r 3 


(4c) 


Eqs. (4) must include the contribution of all dynamic components to 
the motion of the mass center. Since only the array is considered flex- 
ible, all contributions emanate from the array and are included in these 
integral expressions. The additional term "e" which occurs in Likin's 
work is zero under the assumptions here. The Integral terms result in 
higher order terms and will be deleted upon linearizing the following 
equations. The orbital rate is assumed to be very small thus allowing the 
deletion of terms containing time derivatives of the transformation matrix 
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Thus, Eq. (2) muy be expressed, 


F 8 ■■ m S S cV' - 11 J [r' 5 R + ccR']dr 3 - 2q^y J [rgR + otR* ] dr, 
l r„ 


- {q 1 q 1 + q x )P I^R + aR']dr 3 j - ( E ^“(C 1 ^) 


+ [ (G 1 q 1 > <G 1 q 1 ) + 2(0^)^ + ^2 ^ S r l 5 

T . - . . . J 

+ 2[q 2 + (C 1 q 1 )] (r 5 + r 4 a) s + (r* 5 + r^a) ^ 


(5) 


Out: ol! Plane Bending 

Bending in the array Is considered normal to the array. Therefore, 
the third component equation is of interest here. Excluding higher order 
terms, the desired linearized equation becomes upon extracting the third 
component equation from Eq. (5), 


F 3 “ “ (r 3 ^11^10 + C 12^10 + C 13^10^ 


+ r 4 [C 21^10 + C 22^10 + C 23^10 ]) + (r 50 + r 4 a) ( 


( 6 ) 


It should be noted again that the transformation matrix component's 

C 1 will contain a steady contribution and also a contribution which is a 
lj 

function of the vehicle angular velocity disturbance quantities. The 
latter will lead to higher order terms which may be deleted upon expansion 
of the transformation expression thus preserving the- linear nature of the 
equation of motion. 


n 
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The dependent variables of Eq. (6) may be expressed in frequency 
dependent form which is useful in analyzing array frequencies as corre- 
lated to those of the attitude control system. For a linearized system 
it may be assumed that the motion of the craft and array will be harmonic. 
The solutions may be expressed in complex form as a frequency dependent 
magnitude or summation of mode magnitudes times a time dependent exponen- 
tial function. 

The following solutions shall, therefore, be assumed. 

<4> - [<4»>ie Xt < 7 ‘> 

< r 50> - j[4><«hV r 3 ,elt <7b) 

n ' Xt 

(a) = Z [A(X)].lMr,)e AC (7c) 

j=l J 3 


<f“) - 


(7d) 


In the model chosen for this Investigation, no external forcing is 
assumed to exist. Therefore, the term E a consists of the elastic restor- 
ing and damping forces produced by the array. Such forces are functions 
of the deformation variables and have generally the same form as those 
quantities found on the right hand side of Eq. (6). Thus, further speci- 
fication of these terms will produce no additional information for our 
purposes. 

It is possible to derive the desired set of equations by integrating 
Eq. (6) over the array width, then multiplying by the bending mode shape 
^i^ r 3^ anc * ^Ategrat'ing over the array length. The resulting relation is 
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u oat of equations, one for each mode considered, upon which may be 
applied the orthogonality principle r£ normal modes (Ref. 10). The set, 
therefore, produces relations which are uncoupled in the bending modes. 
Thus, Integrating Eq. (6) over r^ yields, 


F 3 dr 4 " A 


|^ m l ^-^(CutQioCX)] + C 12'^10^^ + C 13^10^^) 


+ ^yR5o ( ^)3 j <l> j (r3 ) ] + m 2 [c^tQ^CX)] + C^lQ^CX)] 


1 

+ 4 3 [qJ 0 (X)] + X E^X)!^)] 


( 8 ) 


where the terms m^ and m^ are defined in Appendix B. 

Multiplying Eq. (8) by ^(r-j) and integrating over r^ gives the 
following set of relations. 

1 

F i = * [* ( m l.f 1 [R S0 a)] j b ij + 

+ Vi (- ( c Wio™ ] + c i 2 < W] + 

+ m 2 a i [ C 2l^io^^ + ^22 ^10^ ^ + C 23^10^^ )J 


(9) 


The integral expressions a^, g i# b^ , and c_^_. are found in Appendix B- 
It should be noted that 


j 


h < r 3 ),j, j (r 3 )dr 3 = 1 


0 for i ± j 

1 for i = j 
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Array Torsion 

An additional set of n' equations will now be derived for the elastic 
torsional motion of the array. Suboequently , the six additional relations 
required for analysis will be derived from rigid body considerations. The 
entire set of equations represents the equations of motion for the space- 
craft and contain those vehicle parameters critical to the dynamic behavior 
of the craft. 

The rotational motion of the array element may be examined by con- 
sidering the Euler equation. 

T s - i 8 (10) 


Expanding this relation in terms of the variables previously presented 
and defining the ’'small" element rotation as, 

f S <W 1 




3r, 


-y ~y 

a l + aa 2 


( 11 ) 


allows the expression for Eq. (10) as. 


T S = jj[ S [C 1 ^ + & S ] + [I 3 (C^) + (C 1 q 1 )I S + I S q 2 


+ q 2 X S - (t S C 1 q 1 + I S q 2 )~)e S + [-'[^(C 1 ^ - q^q'j.) ]~ 
+ I 8 (C 1 ^ ~ q 2 C 1 q 1 )’“- - (C 1 q 1 + q,,) *'[I S (C 3 ^ + q^T 
+ (G 1 q 1 + q 2 )V(G 1 q 1 + q 2 )~] 3® - 1%^ 


N 
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+ t (C 1 q 1 ) 4 q 2 ] I B [(0^^ + q 2 ]J (12) 

s s 

The term 1 la the moment of inertia matrix of the element and 3 is the 

rotation which is assumed to be small. < 

Extracting the component equation about the a^ axis and deleting 

higher order term? yields, 





T s ,_1 *1 

I 22^ C 21 q 10 


G 1 a 2 
C 22 q 10 


1 *3 
C 23 q 10 


a)] 


(13) 


For the liuear equation, the solutions may again be assumed to be of the 
form, 

(q$[ Q ) = X[Q^ 0 (X)]e Xt (14a) 

n 1 

(60 « X 2 MA(X)],iMrJe Xt (14b) 

1=1 113 

(T°) - T| e Xt (14c) 


which yields upon substitution into Eq. (13), 


% = [XI S 22 (c^tQ^CX)] + C X 2 [Q^ 0 (X)] 
+ C 23 [qJ 0 (X)] + X^[A(X)] j 4. j (r 3 ))j 


(15) 


Eq. (15) may be multiplied by the mode shape ^(r^) and then successively 
integrated over r^ and then r^ to give, 
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T 1 ■ [ X ( C 21 [Q 10 W)I + C 22 [<) 10 U)J * °23 [< 4H d i 


+ y- S [A(\)],e 
j"l 2 • 


which is a set: of n' equations. In order to complete the equations of 
motion, the rigid body equations must now be developed. 


3. VEHICLE EQUATIONS OF MOTION 


The translational equations for the vehicle may be derived from 


Newton's second law. 


F = M{t} T r' 


which may be expressed in the {&} frame as 


F - MC if. 


It has been assumed that no external forces are applied on the 
vehicle or array, only pure torques. Therefore, both sides of Eq. (18) 
are identically zero. There will be no translational motion of the mass 
center of the vehicle. 

The rotational equation for the vehicle may be derived from the 
Euler relation, 


S _ dH 
" dt 


Expanded, Eq. (19) gives the matrix relation 
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1' - + | j (r 3 + r^)(m s )(r 5 + r 6 >dr 4 dr 3 


+ ^ S 3 S dr^dr 3 


Expressing Eq. (20) in component form yields the relations, 


T 1 " X A + \ \ raSr 3 (if 50 )dr 4 dr - 


■ r ,2 i'3(r c .+r fi ) , \ 

+ ] J m r 4 r 3 cidt 4 dt3 + J J I u — 2 [ ^ J <V r 3 


* | \ 


s ..3 . , 

m r /, r 50 dt 4 dr 3 


* / 

+ m S (r^) 2 adr^dr 3 + I^adr^dr, 

/ J J ~ 


T 3 “ V4 


Solutions may be substituted analogous to those presented in previous 
sections and the integrals evaluated to give, 


C, - flnXCQi 


J«a>] + x 2 f ? [R^a)].^ + mj) 
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1 \ 
♦ + n jjj 


(24) 


I 2 XtQj 0 (X)] + \ 2 [t, [K^ 0 (X)l J P j 

i 1 


+ £ l'A(X)L(t + u .) 
1-1 J 3 3 - 


(25) 


T 3 - (l 3 X[Q 3 0 (X)]) (26) 

These equations along with those derived in previous sections con- 
stitute the necessary set of equations for description of the appendage 
and craft motion. 


IV. PARAMETERS OF IMPORTANCE 


I. PARAMETER CLASSIFICATION 

The craft's equations of motion may now be collected and are expressed 


as ; 

One set of n array bending equations (i = 1 to n) , 


F l- »UWC3 3 + 1> 

•V’ x 

1 

+ X(pm 1 C 33 R' L + m 2 ) t [A(X) ] ,G^. + m^ f C 2i^io^^ 
j“l ' 

+ C 2Z^10^^ + C 23^10^^) “ m l S i ( C 11^10^^ 
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+ <£,[q 2 10 CX>] + 


(27) 


One Bet of n' array torsional equations (1 ■ 1 to n 1 ), 

T 1 " [ X (°iltlio (W3 + U 22 !< llO <X)1 + h 


+ x 2 >: [A(\>] ,c 

j»i 3 


(28) 


One act of three component vehicle rotational equations (p ■ 1 to 3) , 


T p “ V [Q 10 (X)1 + ^ ^ [^ 0 (X)]j[( gj + 

u j“X 


lp 


n 




+ (e 3 + W] 


(29) 


Examination of equations (27) through (29) show that the following 
classification of parameters may be made. 

a. ) Mass dependent parameters m^, m^, U, and various inertia terms 

and Integrals containing such terms . 

b. ) Geometry dependent parameters such as the integrals R and R' , 

and the array length and width. 

c. ) Damping factors and natural mode frequencies. Frequency para- 

meter X, 

d. ) Mode shape dependent parameters and integrals such as a^, 

^ij * c ij ’ e ij * etc * defined in Appendix B. 
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Of thcHc parameter!} the mans and geometry dependent and possibly 
natural frequencies mny be computed prior to the flight. The terms to 
be measured inflight are then damping frequencies, and corresponding mode 
shapes . 

V. DEVELOPMENT OF MEASUREMENT SYSTEM 
1. SYSTEM INSTRUMENTATION 

It is now necessary to qualitatively establish those methods by 
which the inflight measurement may be accomplished. Three methods of 
measurement are Immediately apparent; direct video recording of appendage 
dynamics, strain gages, and mechanical or piezoelectric accelerometers. 
Methods have been demonstrated which allow extraction of mode shapes and 
frequencies for the simple bending of beams using n video recording of 
the vibration. However, the complexity of the onboard system, and com- 
plexity introduced by the bending-torsional combined vibration eliminate 
this technique from further consideration. 

The interpretation of data for the multiple bending and torsional 
modes favors the use of the accelerometer type of instrument as opposed 
to the strain gage. Mode frequencies and array displacements may be found 
from the direct output of the accelerometers by passing the signal through 
a fi*cquency filtering device. 

Mode shape determination will require subsequent analysis of the 
apparent node conditions and amplitudes at each accelerometer station. 

It is essential that a strict time history be established for accurate 
measurement of phase differences between the motions at all accelerometer 
locations, in addition to securing general amplitude and frequency data. 
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Thin will assure appropriate data for the subseiuenl; determination of mode 
shapes . 

2 • BENDING TORSIONU, MODE CONFIGURATIONS 

A preliminary knowledge of the basic mode shape configurations is 
essential in the determination of accelerometer locations. Although pre~ 
else specification of these modes is Impossible, and in fact what we 
desire to measure, preliminary estimates as to their form may promote 
Increased efficiency of the system and aid in the design and positioning 
of each accelerometer. Typical first four bending mode shapes are shown 
in Fig. 3, and the first four torsional mode shapes in Fig. A. 

3. SYSTEM CONFIGURATION 

Although it is unclear at present how many modes will have frequencies 
within the sensitivity range of the attitude control system (which would 
permit an Immediate truncation thus deleting all other modes from further 
consideration) , it may be assumed that the system of accelerometers will 
have the general configuration as shown in Fig. 5, on the array surface 
itself. Contrary to the simplified configuration used in the determina- 
tion of the important parameters, most arrays consist of a supporting 
structure (i.e. pantograph, etc.) plus the array proper. It is, therefore, 
believed that monitoring should include data from the array and also the 
supporting structure for proper analysis of the dynamic characteristics 
of the appendage. The spacing and number of such devices will be depen- 
dent oh tlie number of modes retained after truncation. It is necessary to 
establish an adequate representation of each of the retained modes from 
data received from the system. This dictates certain restraints on 
. ,:\.uele rometer positioning. If possible, accelerometers should not be 
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positioned on any nodal line of the retained modes of either the array or 
the support. This decreases their utility by restricting the number of 
modeB from which data is obtained. The positioning must take into account 
preliminary estimates of node locations in the complex bending and tor- 
sional mode system. Fig. 5 suggests a primary series of accelerometers 
on each side of the elastic-torsional axis and a secondary series which 
along with data from the primary series may allow extraction of the com- 
plex mode locations for the array. Those distances s^ and p,^ are, as 
previously indicated, functions of the specific modes retained. Instrument 
locations on the supporting structure will bii determined by requirements 
similar to those for the array. 

In considering the number of accelerometer positions to be used in 
the primary series for the array a reasonable representation of the highest 
harmonic mode shape may be acquired from data given by 2(N 4-1) or 3(N + 1) 
accelerometer locations depending on the accuracy desired. Here N is the 
value of the highest harmonic mode retained. The secondary series which 
is to be used as a complimentary system to the primary one is estimated 
to require only about fifty percent of this number of positions. 

4. ACCELEROMETER SAMPLING RATE 

The sampling rate Is a function of the frequency of the highest mode 
and should have a value for each accelerometer of about 4f where f is the 
frequency of the mode in cycles per second. The total data rate for the 
array system is then of the order, 

DATA RATE » 8f(N + 1) to 12f(N + 1) 

As an example, let' 9 assume that the fourth harmonic is the highest 
mode to be retained, with an estimated frequency of about 100 H z> 


This would menn that the data rate from the system would need to be, 

DATA RATE “ 4000 to 6000 bits per sec. 

It would be advantageous to reduce this value. However, a detailed 
analysis of expected mode shapes and frequencies would be required to 
properly effect such a reduction or even investigate the feasibility of 
such a reduction. Such an investigation will be included in a later study, 

5. INFLIGHT MANEUVERS 

It is highly desirable to minimize any special inflight maneuvers or 
deviations from the normal flight sequence. The forcing of the array 
system necessary to produce vibrations in the required modes may be 
accomplished by a properly timed firing sequence of the attitude control 
thrusters. This would consist of an Impulse torque applied about a given 
control axis with a subsequent correction torque applied at a specified 
time which would be determined by a modal frequency analysis. Depending 
upon the orientation of the thrusters , the required torques may be applied 
by any appropriate combination of thrusters. No other maneuver is neces- 
sary. 

For data interpretations, the accelerometer system allows reasonably 
straightforward analysis of output. The desired mode frequencies may be 
extracted directly by analysis of the signal from each accelerometer loca- 
tion and the determination of those frequencies to which the given 
accelerometer is responding. The mode shapes will require geometrical 
analysis of the motion at each accelerometer position and the phase 
relationship between these motions. General guidance in this effort will 
be extended by the numerical results of a computer simulation not yet 
completed. 
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VI. CONCLUSIONS 

It Is possible to state several interim conclusions concerning the 
various aspects of the flexibility-control Interaction phenomenon. 
Generally! these interim conclusions represent the nature of important 
parameters, measurement techniques, processing, and data retrieval and 
have not been in all cases rigorously substantiated. However, this pre- 
liminary work has served to Indicate more clearly those areas of investi- 
gation which will produce necessary data for complete design of the 
inflight system which is to be accomplished in follow-on work. It is in 
this context that the following conclusions are presented. 

a. ) Any experiment of this type is very complex and data is not 

unique. 

b. ) The important structural parameters to the flexibility-control 

interaction problem are appendage and vehicle mass and inertia 
characteristics, geometrical properties, vibration damping and 
frequency values, forcing frequencies, and bending and torsional 
mode shapes . 

c. ) Those parameters which must be measured inflight are the damp- 

ing modal frequencies and respective mode shapes. 

d. ) The nature of inflight systems and analysis of system output 

favors the use of accelerometers for the determination of 
frequencies and mode shapes. 

e. ) For the structure flexibility properties, the number of 

accelerometers necessary to create a useful measurement system 
is dictated by the degree of the highest harmonic mode retained 
in the analysis after truncation. It is felt that the number 
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should be oE the order 2(N + 1) to 3(N + 1) for the primary 
series of accelerometers and about half of that value for the 
secondary series. Such speculation requires further sub- 
stantiation using a computer analysis of the combined bend- 
ing-torsion node locations. 

f. ) The data rate for each accelerometer position in a function 

of the frequency of the highest harmonic mode retained. It 
is felt for accurate representation of the modes, a data rate 
of 4f is required where f is the frequency of the highest mode 
in cycles per second. This would indicate a system capability 
requirement of at least 12f(N + 1) data bits per second. It 
is felt that this value can be made feasible in most cases of 
solar array vibration. However, this will impose certain 
design restrictions on the array to be used in conjunction with 
the system design. 

g. ) Special inflight maneuvers will be a sequenced activation of 

the attitude control thrusters for excitation of the vibratory 
modes of the array, both torsional and bending. 

h. ) Data interpretation will consist of a frequency filtering of 

accelerometer output and a geometrical and phase analysis of 
such output for mode shape determination. Specific schemes 
for such interpretation will need to be presented at a later 
time. 

VII. RECOMMENDATIONS 

It is recommended that the following tasks be carried out in a follow 


on study: 
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a. ) A comparison of important parameters for the continuous roll 

out versus the flexibly connected rigid fold out array should 
be performed. Experimental aspects will be fitted to either 
type to the extent possible. 

b. ) The spacecraft attitude control system must be modelled to 

study interaction properties. This will be integrated into 
previous work to develop flexibility /control Interaction 
phenomenon . 

c. ) Computer simulations will be carried out to establish a 

correlation matrix for the important parameters to be 
measured. Techniques for independent measurement will then 
be considered. 

d. ) An analysis of dynamic response characteristics for combina- 

tions of antisymmetric and symmetric bending and torsional 
modes will be developed. 

e. ) Candidate measurement devices and system options will be con- 

sidered to obtain a "best" set up for collecting data. Posi- 
tion and implementation requirements will be included. 

Criteria for evaluation of expected data will also be estab- 
lished. 

Several other tasks may also be considered which include the effects 
of mass discontinuities, a nondimensional study of the equations of motion 
to obtain a general interpretation of results through the use of the laws 
of similitude, and interface and system integration effects on vehicle 
power, structure, TT&C, and propulsion. 
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NEW TECHNOLOGY STATEMENT 


There are no reportable items of new technology. 
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APPENDIX A 


POSITION VECTORS 


The following vectors are necessary to describe the position 
arbitrary element on the solar array. 
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Note: The vector r 2 contains no steady components as {a} and {b} 

assumed to coincide in the undeformed configuration. 
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The angular velocity of the center of «nss with respect to inertial 
space is, 
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while the rotation of the array to maintain proper solar orientation gives, 



(A. 8) 
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One additional matrix form Is useful in expressing matrix equations. 



(A. 9) 


This form is required in expressing cross products such as, 


A 1 x A 2 = ^^ A l x ^'®^ Ta 2 ^ te^A^ 


(A. 10) 
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APPENDIX B 


INTEGRAL EXPRESSIONS 


The following integral expressions occur in the derivation of the 
equations for array bending and torsion, and vehicle rigid body motion. 
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